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ABSTRACT 


This  two-part  series  of  papers  is  concerned  with  the  response  and  various  instabilities 
which  govern  the  behavior  of  circular  cylindrical  shells  under  pure  bending.  Part  I 
describes  the  experimental  part  of  the  study  and  Part  II  presents  the  numerical  sim¬ 
ulation  of  the  various  phenomena  observed  experimentally.  Experiments  were  con¬ 
ducted  on  long  aluminum  6061-T6  shells  with  eleven  different  diameter- to-thickness 
ratios  ranging  from  60.5  to  19.5.  For  such  geometries,  the  structural  response  and  in¬ 
herent  instabilities  are  strongly  influenced  by  the  plastic  characteristics  of  the  material. 

Thinner  shells  were  found  to  develop  short  wavelength  periodic  ripples  on  the 
compressed  side  of  the  shell.  The  shells  buckled  locally  and  collapsed  soon  after  the 
appearance  of  the  ripples.  Thicker  shells  were  found  to  exhibit  a  limit  load  instability 
as  a  direct  consequence  of  the  ovalization  induced  to  the  shell  cross  section  by  bending. 
Following  the  limit  load,  the  ovalization  was  found  to  localize,  leading  to  the  eventual 
collapse  of  the  shells.  For  shells  with  intermediate  Djt  values,  short  wavelength  ripples 
developed  at  the  same  time  that  localization  of  ovalization  was  recorded.  The  shells 
buckled  locally  and  catastrophically  following  the  development  of  a  limit  load. 

INTRODUCTION 

This  study  is  concerned  with  the  problem  of  bending  of  long,  circular  cylindrical  shells. 
Of  particular  interest  is  the  response  of  such  shells,  bent  into  the  plastic  range  of  the 
material,  and  the  various  instabilities  which  limit  the  extent  to  which  the  shells  can  be 
loaded  or  deformed. 

It  is  well  known  that  bending  of  relatively  thin-walled  circular  tubes  induces 
ovalization  to  the  tube  cross  section  (Brazier  effect  (1927)).  The  growth  of  ovalization 
causes  a  progressive  reduction  in  the  bending  rigidity  of  the  shell.  Eventually,  a  maxi¬ 
mum  value  of  moment  is  reached.  Further  bending  occurs  at  a  dropping  moment.  For 
linearly  elastic  shells,  when  the  maximum  moment  is  reached  the  shell  has  undergone 
an  ovalization  which  reduces  the  diameter  of  the  cross  section  in  the  plane  of  bending 
by  approximately  2/9ths.  In  practice,  for  thin  shells,  this  limit  load  instability  is  of¬ 
ten  preceded  by  shell  bifurcation-type  instabilities  characterized  by  circumferential  and 
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axial  waves  (see  Stephens  et  al.  (1975),  Fabian  (1977),  Axelrad  (1980)). 

In  the  case  of  thicker  shells  {D/t  <  100,  say),  the  response,  as  well  as  the  insta¬ 
bilities,  are  strongly  influenced  by  interaction  between  the  induced  ovalization  and  the 
plastic  characteristics  of  the  material.  Ades  (1957),  calculated  the  nonlinear  moment- 
curvature  response  of  long  elastic-plastic  tubes  undergoing  uniform  ovalization.  He 
assumed  the  cross  section  to  deform  into  an  ellipse  and  used  the  deformation  the¬ 
ory  of  plasticity  to  model  the  nonlinear  material  response.  Gellin  (1980),  developed 
a  more  accurate  solution  of  the  problem  using  improved  kinematics  and  the  defor¬ 
mation  theory  of  plasticity.  Gellin ’s  formulation  and  solution  procedure  were  further 
refined  by  Shaw  and  Kyriakides  (1985)  by  adopting  a  more  complete  set  of  extensional 
kinematics  and  incremental  plasticity  models  (see  also  Fabian,  1981,  Bushnell,  1981). 

These  solutions  demonstrated  that,  although  the  ovalization  induced  to  the  shell 
cross  section  is  relatively  small  for  shells  bent  into  the  plastic  range,  it  still  leads  to 
the  development  of  a  limit  moment  in  the  response.  This  is  demonstrated  in  Fig.  1, 
where  three  moment-curvature  responses  for  an  aluminum  shell  with  a  diameter-to- 
thickness  ratio  {Djt)  of  35.7  are  compared.  One  of  the  responses  shown  was  obtained 
by  simple  integration  of  the  stresses  through  the  cross  section  (Saint- Venant  bending); 
as  a  result  the  moment  is  monotonically  increasing.  The  second  response  was  calculated 
by  a  Brazier-type  analysis,  which  includes  the  effect  of  the  ovalized  cross  section.  It 
has  a  somewhat  lower  moment  but,  more  importantly,  the  ovalization  results  in  the 
development  of  a  limit  moment. 

Experiments  involving  pure  bending  of  inelastic  shells  have  been  carried  out  by 
Moore  and  Clark  (1952),  Jirsa  et  al.  (1972),  Sherman  (1976),  Tuggu  and  Shroeder 
(1979),  Reddy  (1979),  Kyriakides  and  Shaw  (1987),  and  others.  Johns  et  al.  (1979), 
Corona  and  Kyriakides  (1988)  and  Ju  and  Kyriakides  (1990)  carried  out  experiments 
involving  combined  bending  and  external  pressure. 

The  experiments  demonstrated  that,  in  addition  to  the  limit  load  instabilities, 
shell  bending  is  limited  by  various  shell  buckling  modes.  In  addition,  evidence  of  dis¬ 
tinct  localization  of  deformation  was  observed  in  many  shells  (see  Fig.  6  in  Corona  and 
Kyriakides,  1988).  Figure  1  includes  the  moment-curvature  response  obtained  from  an 
experiment  on  an  aluminum  606TT6  shell  with  Djt  =  35.7.  The  shell  was  found  to 
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ovalize  uniformly  for  most  of  the  bending  history.  At  the  curvature  indicated  by  (|) 
short  wavelength  ripples  were  observed  on  the  compressed  side  of  the  shell.  Soon  after 
the  ripples  appeared,  the  ovalization  started  to  localize  in  a  region  approximately  8 
shell  diameters  long,  a  hmit  load  developed  on  the  response  (identified  by  “A”)  and 
the  shell  buckled  locally  at  a  dropping  moment.  A  buckled  test  specimen  is  shown  in 
Fig.  2a.  As  the  hmit  load  was  approached,  it  was  observed  that  long  wavelength  initial 
imperfections  in  the  shell  tended  to  be  amplified.  This  resulted  in  some  nonuniform 
growth  of  ovaUzation.  An  example  of  this  behavior  is  shown  in  Fig.  3  where  axial 
scans  of  the  change  in  diameter  in  the  shell  (ADfDo),  in  the  plane  of  bending,  taken 
at  different  values  of  curvature  are  shown.  The  shell  cross  section  ovahzed  in  a  rea¬ 
sonably  uniform  manner  up  to  a  curvature  of  k/ki  =  0.80.  At  higher  curvatures,  long 
wavelength  axial  waves  are  seen  to  develop.  The  shell  failed  catastrophically  by  local 
buckling  in  one  of  the  troughs  of  ovalization. 

Clearly  the  behavior  described  above  is  different  and  more  complicated  than  what 
has  been  modeled  in  the  references  mentioned  above.  The  objective  of  this  study  was 
to  identify  and  understand  the  role  of  the  major  instabilities  which  affect  the  response 
of  long  elastic-plastic  shells  under  pure  bending.  This  was  done  through  a  combination 
of  experiment  and  analysis.  A  systematic  set  of  experiments  were  conducted  using 
aluminum  shells  with  a  broad  range  of  diameter- to-thickness  ratios.  The  experimental 
facilities  used,  procedures  followed,  the  results  from  the  experiments  and  a  discussion 
of  the  observations  are  presented  in  part  I.  The  major  phenomena  observed  in  the 
experiments  are  simulated  numerically  in  part  11.  The  numerical  predictions  are  used 
to  explore  and  explain  the  behavior  of  the  structures  tested. 

EXPERIMENTS 
a.  Test  Facilities 

The  bending  experiments  were  conducted  in  a  test  facility  consisting  of  a  pure  bending 
device,  shown  in  Fig.  4,  various  transducers  for  recording  geometric  changes  in  the 
shell  and  a  computer- b2ised  data  acquisition  system.  The  bending  device  was  designed 
as  a  four-point  bending  machine  capable  of  applying  bending  and  reverse  bending  (see 
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Kyriakides  and  Shaw  -  1987).  It  consists  of  two  double  strand  sprockets  (pitch  diameter 
9.56  in  (243  mm))  mounted  with  pillow  blocks  on  two  parallel  beams  as  shown  in  Fig  5. 
The  beams  are  70  in  (1.78  m)  long  and  the  distance  between  the  sprockets  can  be 
varied.  In  the  experiments  described  here  the  sprockets  were  located  49.5  in  (1.26  m) 
apart.  Heavy  chains  run  around  the  sprockets  and  are  connected  to  hydrauhc  cyhnders 
and  load  cells  forming  a  closed  loop.  For  monotonic  bending,  the  bending  curvature 
capacity  of  the  machine  can  be  doubled  by  removing  one  of  the  cylinders  and  extending 
the  second  one  to  its  full  stroke  of  10  in  (0.25  m). 

Each  test  specimen  is  fitted  with  close-fitting  solid  extension  rods,  as  shown  in 
Fig.  5.  The  rods  extend  more  than  3  diameters  into  the  free  ends  of  the  shell  and 
are  carefully  chamfered  at  the  ends  to  reduce  the  effect  of  the  discontinuity.  The  test 
specimen  assembly  engages  the  sprockets  through  four  rollers,  as  shown  in  the  figure. 

Bending  is  achieved  by  contracting  one  of  the  hydraulic  cylinders,  causing  rotation 
of  the  sprockets.  During  bending,  the  rolling  contact  between  the  end  rods  and  the  test 
facility  allows  free  axial  movement  of  the  test  specimen.  While  the  deformation  of  the 
test  specimen  is  uniform  along  the  length,  the  bending  curvature  is  directly  proportional 
to  the  rotation  of  the  two  sprockets.  The  rotation  of  each  sprocket  is  monitored  by  a 
linear  voltage  differential  transformer  (LVDT)  as  follows;  A  flexible,  inextensional, 
thin  cable  runs  over  the  circular  hub  of  the  sprocket.  The  cable  is  spring-loaded  and  is 
connected  to  the  core  of  the  LVDT.  Thus,  rotation  causes  linear  movement  of  the  core 
and  a  proportional  signal  from  the  LVDT. 

The  majority  of  the  test  specimens  tested  had  an  effective  length  of  30  in  (0.76  m). 
With  the  test  facility  configuration  used  such  specimens  could  be  bent  to  a  curvature 
of  0.070  in~^  (2.74  m“*).  Larger  curvatures  can  be  achieved  by  reducing  the  length  of 
the  shell. 

The  bending  moment  in  the  test  specimen  is  directly  proportional  to  the  tension 
in  the  chain.  This  is  monitored  by  calibrated  load  cells.  The  capacity  of  the  bending 
device  is  15,000  lb-in  (1,700  N-m).  The  device  was  designed  to  be  stiff  relative  to  the 
test  specimens.  For  the  experiments  reported  here,  the  energy  stored  in  the  device  was 
less  than  4.5%  of  the  energy  stored  in  the  strongest  of  the  test  specimens  (approximately 
half  of  this  is  stored  in  the  hydraulic  system).  A  stiff  machine  is  essential  for  testing 
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structures  beyond  limit  loads. 

The  ovalization  of  the  cross  section  of  the  tube  is  known  to  play  an  important 
role  in  the  stability  of  the  structure.  Thus,  several  transducers  for  monitoring  small 
changes  in  the  geometry  of  the  shell  during  bending  were  developed.  Figure  6  shows 
a  lightweight  transducer  which  is  mounted  on  the  shell  and  used  to  monitor  changes 
in  the  diameter  in  the  plane  of  bending.  It  consists  of  a  rigid  frame  which  supports 
a  sliding  block.  The  instrument  is  lightly  spring  loaded,  as  shown  in  the  figure,  and 
contacts  the  shell  through  knife  edges.  The  relative  displacement  between  the  two  knife 
edges  is  monitored  by  a  miniature  LVDT.  In  some  experiments,  two  such  transducers, 
placed  approximately  ten  shell  diameters  apart,  were  used  in  order  to  establish  the 
onset  of  localized  deformation  in  the  shell. 

A  computer-based  data  acquisition  system  was  used  to  monitor  the  variables 
measured  during  the  experiment.  It  consists  of  several  signal  conditioning  units,  an 
analog-to-digital  (A/D)  converter,  a  desktop  computer  and  peripherals  such  as  digital 
plotters  and  an  electronic  display  unit.  The  system  enables  real  time  display  of  some 
of  the  variables  and  provides  storage  of  the  data  for  later  analysis. 

During  a  typical  monotonic  bending  experiment  on  a  shell,  the  variables  moni¬ 
tored  were  the  moment  (M),  the  curvature  (k)  and  the  ovalization  at  one  or  two  discrete 
points  along  the  length  of  the  shell  (AD,).  The  data  acquisition  system  stores  a  set  of 
{M,k,  AA}j  whenever  any  one  of  the  variables  changes  by  a  prescribed  small  value. 
The  number  of  sets  of  data  collected  from  each  experiment  varied  from  approximately 
125,  for  the  thinner  shells  to  450  for  the  thicker  shells.  The  uncertainties  in  the  recorded 
data  were  as  follows:  moment  ~  0.8%;  curvature  ~  0.6%;  ovalization  ~  0.1%. 

In  at  least  one  experiment  for  each  shell  Djt  considered,  a  second  transducer, 
shown  in  Fig.  7,  was  used  to  periodically  scan  the  ovalization  along  the  length  of  the 
shell.  In  this  case,  contact  with  the  shell  is  made  through  rollers.  The  transducer  is 
moved  along  the  length  of  the  shell  manually.  Its  axial  position  is  established  by  an 
encoder  system  which  works  as  follows:  a  thin  flexible  tape,  with  regularly  spaced  black- 
white  markings,  is  bonded  on  the  tension  side  of  the  shell  as  shown  in  the  figure.  An 
emitter/receiver  photodiode,  mounted  on  the  transducer,  produces  an  electrical  pulse 
whenever  it  encounters  a  white  marking.  The  pulse,  suitably  processed,  triggers  the 
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A/D  converter  and  the  value  of  ovalization  is  recorded  in  the  computer.  In  the  majority 
of  the  experiments  presented  here,  the  encoder  had  a  resolution  of  five  readings  per  inch. 

In  the  case  of  thinner  shells,  short  wave  length  ripples  develop  on  the  compression 
side  of  the  bent  shell.  For  the  aluminum  shells  used  in  the  experiments,  with  typical 
diameters  between  1.0  and  1.5  in  (25-38  mm),  the  amplitude  of  these  ripples  was  of  the 
order  of  0.001  -  0.005  in  (0.025  -  0.125  mm).  Thus,  a  more  sensitive  instrument  with 
higher  axial  resolution,  shown  in  Fig.  8,  was  developed  for  measuring  the  geometry  of 
these  ripples.  It  consists  of  an  axial  guide,  with  a  fine  lead  screw,  which  supports  a 
vertical  traversing  system.  This  is  lightly  spring-loaded  and  is  in  contact  with  the  top 
and  bottom  of  the  shell  through  small  diameter  rollers.  The  guide  is  4  in  (102  mm)  long 
and  is  mounted  on  the  bent  shell  as  shown  in  the  figure.  The  vertical  frame  is  moved 
axially  by  turning  the  lead  screw.  The  two  traversing  beams  are  free  to  move  vertically 
to  accommodate  changes  in  the  geometry  of  the  shell.  Changes  in  the  distance  they 
are  apart  are  monitored  by  a  proximity,  non-contact,  displacement  transducer  with  a 
sensitivity  of  125  mV/10“^  in  (4.92  V/mm).  The  axial  position  of  the  transducer  is 
monitored  by  an  encoder  system  which  consists  of  an  emitter/receiver  photodiode  and 
a  circular  disc  with  six  equally- spaced  narrow  slits.  The  photodiode  produces  a  pulse 
whenever  the  light  beam  is  obstructed.  The  signal  from  the  encoder  is  used  to  trigger 
the  A/D  converter  and  the  computer  stores  the  output  of  the  proximity  transducer. 
The  axial  resolution  of  the  system  provided  for  192  readings  per  inch. 

b.  Experimental  Procedure 

The  experiments  were  carried  out  on  commercially  available  aluminum  6061 -T6  drawn 
tubes.  The  test  specimens  had  nominal  diameters  ranging  from  1.0  in  (25.4  mm)  to  1.5 
in  (38.1  mm).  Tubes  with  eleven  different  diameter-to-thickness  ratios,  ranging  from 
60.5  to  19.5,  were  tested.  The  lengths  (effective)  of  the  test  specimens  were  30  in  (0.76 
m);  thus,  their  length-to-diameter  ratios  ranged  between  18  to  30.  Average  values  of 
the  diameters  and  wall  thicknesses  of  the  test  specimens,  as  well  as  their  lengths,  are 
listed  in  Table  1.  In  each  case,  enough  measurements  were  made  to  also  establish  the 
maximum  value  of  initial  ovalization  along  the  length  of  the  tube.  The  maximum  value 


6 


of  ovalization  measured,  Aq  is  defined  as  follows: 


Ao  — 


DmAX- Dr, 


Dr, 


(1) 


'max  d'  Drtun 

The  Vcilues  of  Ao  measured  in  each  shell  are  hsted  in  Table  1 .  Wall  thickness  variations 
measured  at  the  two  ends  of  the  test  specimens  are  also  given  by  their  extreme  value 
Ho  defined  as  follows: 

/  _  #  . 
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The  material  properties  of  each  tube  tested  were  measured  from  uniaxial  tests 
on  axial  specimens  cut  from  the  same  length  of  tube  as  the  bending  specimens.  The 
stress-strain  curves  obtained  were  fitted  with  the  Ramberg-Osgood  fit  given  by 


c 
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The  fit  parameters  calculated  from  each  test  specimen  are  listed  in  Table  1.  In  addition, 
the  0.2%  strain  offset  yield  stress  (ao)  is  listed  for  each  case.  In  a  number  of  the  test 
specimens,  an  additional  experiment  was  conducted  through  which  the  yield  stress  in 
the  circumferential  direction  was  established  (as  in  Kyriakides  and  Yeh,  1988).  The 
variable  S  represents  the  ratio  of  the  yield  stress  in  the  circumferential  to  that  in  the 
axial  direction  (the  yield  stress  in  the  thickness  direction  was  assumed  to  be  the  same 
as  that  in  the  axial  direction). 

A  minimum  of  two  experiments  were  conducted  for  each  Dji  considered.  The 
first  involved  monotonic  bending  of  the  test  specimen  at  a  strain  rate  (maximum)  of 
approximately  1%/min.  Each  specimen  was  loaded  to  failure. 

In  the  second  experiment,  the  bending  was  periodically  interrupted.  The  flow'  to 
the  hydraulic  cylinder  was  shut  and  then  an  ovalization  scan  was  conducted.  Six  to 
twelve  such  scans  were  taken  during  the  bending  history.  During  every  interruption 
of  the  loading,  the  compressed  side  of  the  shell  was  visually  inspected  for  the  possible 
appearance  of  ripples.  When  such  ripples  were  detected,  the  instrument,  shown  in 
Fig.  8,  was  installed  on  the  test  specimen  and  used  to  conduct  a  detailed  examination 
of  the  area  in  which  they  appeared.  Ripples  usually  appeared  in  pockets  of  two  to 
four  randomly  distributed  along  the  length.  Figure  9a  shows  a  picture  of  two  such 
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pockets  of  ripples  for  test  specimen  No.  7.  The  range  of  curvature  during  which  the 
ripples  were  visible  was  rather  limited.  Soon  after  their  appearance,  the  shell  buckled 
catastrophically,  developing  one  distinct  local  buckle.  An  example  of  such  a  buckle  is 
shown  in  Fig.  9b.  Following  such  a  buckle,  the  bending  capacity  of  the  shell  usually 
dropped  by  80  or  90%. 

RESULTS 

The  experimental  results  will  be  presented  in  nondimensional  form  using  the  following 
normalizing  parameters: 

Mo  =  (ToDlt  ,  /Cl  =  t/Dl  .  (4) 

The  curvature  at  which  ripples  were  first  detected  on  the  compressed  side  of  the  shell 
is  identified  by  Kb,  in  table  2,  and  by  the  symbol  (J.)  in  the  figures.  (Note  that  /cj  was 
usually  obtained  from  the  second  bending  experiment  in  which  bending  was  interrupted 
while  the  shell  was  checked  for  ripples.  Thus  the  search  was  discrete  and  as  a  result 
the  values  of  Kb  recorded  are  in  general  somewhat  larger  than  the  actual  values.)  The 
curvature  at  which  the  maximum  moment  was  achieved  (limit  moment)  is  identified  by 
Ki  in  table  2  and  by  the  symbol  (A)  in  the  figures.  The  curvature  at  which  catastrophic 
collapse  occurred  is  identified  as  Kc  and  by  the  symbol  (J.).  For  completeness,  the  max¬ 
imum  Saint-Venant  bending  strains  corresponding  to  each  of  these  critical  curvatures 
(cfc,  £l,  €c)  are  also  included  in  the  table. 

The  analysis  developed  by  Shaw  and  Kyriakides  (1985)  was  used  to  calculate 
the  moment-curvature  and  ovalization-curvature  responses  for  each  of  the  experiments 
conducted.  In  this  analysis,  the  shell  is  treated  as  a  long  tube  whose  cross  section 
can  ovalize  but  is  otherwise  uniform  along  the  axis.  The  analysis  is  thus  capable  of 
predicting  the  limit  load  instability  inherent  in  the  problem.  The  predictions  from  this 
analysis  are  shown  by  a  dcished  line  in  each  of  the  figures. 

The  moment-curvature  and  ovalization-curvature  responses  from  the  test  speci¬ 
men  with  the  largest  value  of  0(1  tested  {D/t  =  60.5)  are  shown  in  Fig.  10.  The 
shell  buckled  catastrophically  at  a  curvature  of  KcJ k\  =  0.93  which  is  equivalent  to  a 
maximum  Saint-Venant  bending  strain  of  0.79%.  Buckling  occurred  at  an  increasing 
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moment  and  took  the  form  of  one  sharp,  local  buckle  characterized  by  a  number  of 
diamond  shapes  on  the  compressed  side  of  the  shell,  similar  to  the  one  shown  in  Fig. 
9b  (the  diamond  shapes  are  reminiscent  of  elastic  buckling  of  shells).  The  ovalization 
measured  along  the  length  of  the  shell,  at  various  values  of  curvature,  are  shown  in  Fig. 
11.  Apart  from  the  boundaries,  the  ovalization  is  seen  to  remain  essentially  uniform 
along  the  length  for  all  values  of  curvature  shown.  For  this  value  oi  D  ft,  no  axial  ripples 
were  detected  for  either  of  the  shells  tested. 

The  moment-curvature  «ind  ovalization-curvature  responses  predicted  numerically 
are  seen  to  be  in  very  good  agreement  with  the  experimental  ones.  The  curvature  at 
which  the  limit  load  was  predicted  to  occur  is  seen  to  be  approximately  90%  higher 
than  the  collapse  curvature. 

The  behavior  of  the  shell  with  Djt  =  52.6  was  very  similar  to  that  described 
above.  The  curvature  and  the  bending  strain  at  collapse  are  listed  in  table  2.  The 
major  results  from  the  shells  with  D/i  values  of  50  and  44  are  shown  in  Figs  12-15. 
The  major  features  of  the  responses  are  essentially  the  same.  As  the  D/i  decreases, 
the  Saint- Venant  bending  strain  and  the  curvature  at  collapse  increase.  In  addition, 
the  difference  between  the  curvature  corresponding  to  the  predicted  limit  moment  and 
that  at  collapse  decreases. 

Axial  ovalization  scans  taken  at  different  values  of  curvature  are  shown  in  Figs  13 
and  15.  The  ovalization  is  seen  to  remain  reasonably  uniform  throughout  the  bending 
history  for  both  shells.  In  both  cases,  pockets  of  short  wave  length  ripples  were  detected 
on  the  compressed  side  of  the  shell  just  prior  to  catastrophic  collapse.  The  location  of 
one  of  the  pockets  is  identified  in  the  ovalization  scans.  A  detailed  axial  scan  of  the 
ripple  pocket  is  shown,  in  each  case,  in  the  inset.  The  average  wavelength  of  the  ripples 
measured,  A/\/^,  is  given  in  table  2. 

Although  the  appearance  of  the  ripples  during  the  bending  history  was  relatively 
short,  in  some  cases  it  was  possible  to  monitor  their  development  by  conducting  local, 
axial  scans  of  the  shell  at  small  intervals  of  curvature.  As  the  curvature  increased  the 
ripples  grew  in  amplitude  but  remained  fixed  at  their  original  position  along  the  length 
of  the  shell.  During  the  course  of  this  study  it  was  found  that  improperly  designed  load 
transfer  mechanisms  at  the  boundaries  can  result  in  premature  buckling  of  the  shell 


9 


with  the  buckle  developing  at  the  ends. 

Figures  16  and  17  show  the  results  from  a  shell  with  D ft  =  35.7.  For  this  case,  and 
for  all  other  shells  tested  with  lower  Djt  values,  the  axial  scans  of  ovalization  at  different 
values  of  curvature  produced  evidence  that,  beyond  a  certain  value  of  curvature,  the 
ovalization  started  growing  in  a  non-uniform  fashion  along  the  length.  For  example,  in 
Fig.  17,  it  can  be  seen  that  the  ovalization  remained  uniform  up  to  a  value  of  /c//ci  = 
0.81.  The  beginning  of  some  non-uniform  growth  is  seen  to  occur  at  /c//ci  =  0.93. 
At  the  last  value  of  curvature  at  which  an  axial  scan  was  taken,  k/ki  =  0.99  (i.e., 
just  prior  to  collapse),  the  ovalization  is  seen  to  grow  significantly  more  in  a  region 
approximately  8  diameters  long.  The  maximum  ovalization  recorded  is  approximately 
twice  the  average  value  over  the  remainder  of  the  shell.  At  curvature  of  Kb/Ki  =  0.96, 
a  pocket  of  short  wavelength  ripples  be.  ■'  visible  in  the  trough  of  the  o\'erly  ovalized 
region.  A  detailed  axial  scan  of  this  region  is  shown  in  the  inset  in  Fig.  17.  Further 
bending  resulted  in  catastrophic  buckling.  The  buckle  was  triggered  by  the  ripples,  and 
had  the  form  of  one  sharp,  local  kink,  similar  to  the  one  shown  in  Fig.  2a. 

The  moment-curvature  and  ovalization-curvature  responses  recorded  for  this  shell 
are  shown  in  Fig.  16  (different  test  specimen  from  that  in  Fig.  17).  The  moment- 
curvature  response  developed  a  limit  moment  at  a  curvature  of  kl/kj  =  0.99.  The 
shell  collapsed  catastrophically  at  Kc/kj  =  1.09.  Due  to  the  evidence  of  some  localized 
deformation,  two  ovalization  transducers,  placed  7.2  shell  diameters  apart,  were  used 
in  this  experiment.  Their  axial  positions,  relative  to  the  position  of  the  final  buckle, 
are  shown  in  the  inset  in  Fig.  16b.  The  measured  ovalizations  at  the  two  points  are 
seen  to  maintain  the  same  value  up  to  a  curvature  of  approximately  k/ki  =  0.9.  At 
higher  values  of  curvature,  the  ovalization  is  seen  to  grow  faster  at  point  B  than  at 
point  A.  The  proximity  of  point  B  to  the  position  of  the  buckle,  and  the  results  from 
the  experiment  in  Fig.  17,  explain  this  behavior. 

We  observe  that  the  limit  load  predicted  by  assuming  the  shell  to  remain  uniform 
along  its  length  is  significantly  higher  than  the  actual  limit  load  recorded  in  the  exper¬ 
iment.  We  thus  conclude  that  the  cause  of  the  two  instabilities  is  probably  different. 
We  point  out  once  more,  that  the  cause  of  failure  was  again  short  wavelength  buckling. 
Quite  clearly,  some  interaction  between  the  observed  localization  in  ovalization  and  the 
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short  wavelength  buckling  mode  does  take  place.  It  would  be  of  interest  to  understand 
this  interaction  and  its  effect  on  the  final  curvature  at  collapse. 

Figure  18  shows  the  structural  response  from  a  shell  with  D jt  =  32.2.  The  main 
characteristics  of  the  results  are  similar  to  those  of  the  shell  with  Djt  =  35.7.  Because 
of  the  thicker  wall,  the  difference  between  the  measured  and  calculated  curvatures  at  the 
corresponding  limit  moments  is  smaller.  Short  wavelength  ripples  were  again  observed 
prior  to  the  limit  load.  Figure  19  shows  four  local  axial  scans  of  the  diameter  of  the  shell 
in  the  plane  of  bending  taken  in  the  curvature  range  o{  k/ki  =  0.934  and  «/ki  =  0.991. 
The  results  illustrate  the  growth  of  a  pocket  of  two  ripples.  The  amplitute  of  the  ripples 
is  seen  to  grow  as  the  curvature  was  increased;  but  their  axial  position  remained  fixed. 

In  the  neighborhood  of  the  limit  load  the  shell  again  developed  some  localized 
ovalization  (see  Fig.  18b).  The  cause  of  the  observed  localization  will  be  explored  in 
the  following  paper.  At  this  stage,  it  suffices  to  say  that,  for  shells  with  intermediate 
values  of  D/t  (i.e.,  approximately  26  <  D/t  <  40),  long  wavelength  imperfections  tend 
to  be  amplified,  to  some  degree,  as  the  collapse  curvature  is  approached.  The  result 
is  localized  growth  of  ovalization,  as  shown  in  Figs  3  and  17.  Similarly,  areas  with 
some  stress  concentration,  such  as  might  be  caused  by  overly  constraining  boundary 
conditions,  can  also  cause  some  non-uniform  growth  of  ovalization. 

A  set  of  results  obtained  from  a  shell  with  Djt  =  25.3  is  shown  in  Figs  20  and  21. 
The  measured  moment-curvature  response  is  shown  in  Fig.  20a.  Following  yielding,  the 
moment  is  seen  to  follow  a  protracted  plateau  (very  similar  to  stress-strain  response) 
where  bending  occurs  at  a  mildly  increasing  moment.  At  a  curvature  of  Kij k\  =  0.95, 
the  moment  reaches  its  maximum  value.  The  corresponding  Saint-Venant  bending 
strain  is  2.03%.  Following  the  maximum,  the  moment  is  seen  to  drop  significantly  with 
curvature.  The  shell  buckled  at  a  curvature  of  =  105. 

The  ovalizations  measured  at  two  points  along  the  shell  are  shown  in  Fig.  20b. 
The  two  measurements  remained  essentially  the  same  up  to  a  curvature  of  approxi¬ 
mately  k/ki  =  0.6.  At  higher  values  of  curvature,  transducer  A,  which  was  later  found 
to  be  located  close  to  the  position  of  the  final  buckle,  started  to  record  progressively 
higher  values  of  ovalization  than  transducer  B.  Beyond  the  limit  moment,  the  ovaliza¬ 
tions  recorded  by  the  two  transducers  diverged  significantly.  At  point  B,  the  ovalization 
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eventually  started  to  decrease,  indicating  local  unloading  of  the  shell,  whereas  at  point 
A,  the  ovalization  started  to  show  signs  of  unstable  growth.  The  trough  of  the  final 
buckle  was  1.2  diameters  from  point  A.  The  buckle  was  characterized  by  less  abrupt 
gradients  than  those  of  the  buckles  in  thinner  shells.  In  addition,  the  onset  of  buckling 
was  more  gradual  and  the  drop  in  the  post-buckling  moment  bearing  capacity  relatively 
smaller.  However,  the  buckled  shell  can  still  be  considered  to  be  structurally  destroyed, 
in  this  case  no  short  wavelength  ripples  were  detected  anywhere  on  the  shell.  This  was 
also  the  case  for  all  shells  tested  with  Djt  <  25.3. 

The  moment-curvature  and  ovalization-curvature  responses,  calculated  by  the 
uniform  tube  analysis,  are  included  in  the  figures.  The  predicted  hmit  moment  is  12% 
higher  than  the  measured  value,  and  the  ovalization-curvature  response  is  generally  in 
good  agreement  with  the  value  measured  at  point  B  up  to  the  limit  moment. 

Ovalization  scans  along  the  length  of  a  second  shell,  taken  at  various  values  of 
curvature  during  the  bending  history,  are  shown  in  Fig.  21.  In  spite  of  the  obvious 
presence  of  some  initial  geometric  imperfections  along  the  shell,  the  ovalization  is  seen 
to  remain  uniform  along  the  length  up  to  a  curvature  of  approximately  ac/«i  =  0.67. 
At  kJ K\  =  0.73,  a  section  in  the  middle  of  the  shell,  initially  approximately  5  diameters 
long,  started  to  experience  greater  ovalization  than  the  rest  of  the  shell.  At  higher  values 
of  curvature,  the  length  of  this  section  grew  until,  at  a  curvature  of  0.89,  it  reached  a 
length  of  approximately  10  diameters.  At  this  curvature,  the  value  of  the  maximum 
ovalization  wa^  approximately  2.7  times  higher  than  that  in  the  section  of  shell  which 
remained  uniform.  The  axial  position  at  which  the  shell  localized  varied  from  one 
specimen  to  another  depending  on  the  position  of  the  most  favorable  imperfections. 

In  the  case  of  test  specimens  which  undergo  significant  localized  deformation 
beyond  the  limit  moment,  the  global  deformation  variable,  >c,  is  inappropriate.  The 
rotations  of  each  end  of  the  shell  and  the  axial  position  of  the  area  of  localization  can 
be  used  to  analyze  the  post  limit  load  response.  However,  this  would  complicate  the 
presentation  of  the  results.  Thus  for  convenience,  we  continue  to  use  k  as  the  measure 
of  deformation  beyond  the  limit  load.  The  reader  is  cautioned  that  this  curvature  is 
erroneous  to  some  degree.  For  the  same  reason,  the  curvature  at  collapse,  presented  in 
table  2,  for  test  specimens  7  to  11,  should  only  be  used  for  qualitative  comparisons. 
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Results  from  a  shell  with  Djt  of  19.5  are  shown  in  Figs  22  and  23.  In  this  case,  the 
limit  moment  occurred  at  kl/ki  =  0.96,  which  corresponds  to  a  bending  strain  of  2.79%. 
It  is  interesting  to  observe  that  the  curvature  at  the  limit  moment  calculated  by  the 
uniform  tube  analysis  is  only  5%  higher  than  this  value.  The  ovalization  measurements, 
presented  in  Fig.  22,  show  the  values  measured  at  points  A  and  B  to  remain  very  close 
up  to  the  attainment  of  the  limit  moment,  beyond  which  they  diverge  significantly. 
As  a  result  of  the  dropping  moment,  point  A,  which  was  6.5  diameters  away  from 
the  trough  of  localization,  experienced  a  significant  decrease  in  the  rate  of  growth  of 
ovalization  and,  eventually,  the  ovalization  started  decreasing  (a  point  further  away 
from  the  buckle  would  experience  reduction  in  AD  earlier). 

Figure  23  shows  the  axial  ovalization  scans  taken  at  different  values  of  curvature. 
The  shell  is  observed  to  remain  uniform  up  to  a  curvature  of  approximately  k/ki  =  0.75. 
At  higher  curvatures,  a  section  approximately  8  diameters  long,  located  in  the  middle 
of  the  shell,  starts  to  ovalize  at  a  faster  rate.  The  rate  of  growth  in  this  region  increases 
significantly  for  curvatures  higher  than  k/ki  =  0.94.  The  length  of  the  localized  region, 
at  the  last  curvature  at  which  a  scan  was  taken,  was  approximately  10  diameters. 
After  significant  deformation  in  the  post  limit  load  regime,  the  shell  collapsed.  The 
post-collapse  moment  bearing  capacity  was  approximately  50%  of  the  limit  moment. 
The  collapsed  section  is  shown  in  Fig.  2b. 

DISCUSSION  OF  THE  RESULTS 

The  experimental  results  presented  demonstrate  that  the  problem  is  rich  in  phenomena 
which  are  dependent  on  the  D/t  of  the  shell.  The  problem  has  an  inherent  limit  load 
instability  which  governed  the  response  of  some  of  the  shells  tested.  In  addition,  short 
wavelength  buckling  modes  were  identified  for  some  shells.  In  some  cases,  localization 
with  characteristic  lengths  of  the  order  of  8  to  10  shell  diameters  was  found  to  interact 
with  short  wavelength  buckling  modes. 

In  order  to  highlight  the  instabilities  governing  the  behavior  observed  and  identify 
the  goals  of  the  analysis,  we  divide  the  shells  into  the  following  broad  categories; 


13 


a.  D/t>40^ 


In  this  category  of  relatively  thinner  shells,  the  prevalent  mode  of  instability  wcis  short 
wavelength  rippling.  The  curvature  (or  bending  strain)  at  which  the  ripples  appeared 
was  significantly  smaller  than  that  corresponding  to  the  limit  load  resulting  from  uni¬ 
form  ovalization.  As  a  result,  the  two  instabilities  did  not  show  any  signs  of  interaction. 
The  shells  were  found  to  ovalize  quite  uniformly  up  to  collapse. 

The  short  wavelength  ripples  appeared  in  small  pockets  randomly  distributed 
along  the  length.  They  appeared  at  an  increasing  moment  shortly  before  the  shell 
collapsed.  The  amplitude  of  the  ripples  grew  in  a  non-uniform  fashion  (localized).  The 
shells  failed  catastrophically  when  a  second  instability,  triggered  in  one  of  the  ripple 
pockets,  occurred.  This  happened  at  a  slightly  higher  value  of  curvature  than  kj.  The 
collapse  mode  was  characterized  by  diamond  shapes,  common  in  elastic  shell  buckling. 
The  post  buckling  strength  of  the  shells  was  found  to  be  insignificantly  small  to  be  of 
structural  use. 

We  propose  that  the  observed  rippling  qualitatively  corresponds  to  the  axisymmet- 
ric  buckling  mode  common  in  plastic  buckling  of  axially  loaded  cylindrical  shells  (see 
Lee  (1962),  Batterman  (1965)  and  Batterman  and  Lee  (1966)).  For  thinner  shells,  this 
is  followed  by  a  non-axisymmetric  buckling  mode  which  leads  to  catastrophic  collapse 
(see  Gellin  (1979)  and  Tvergaard  (1983)).  The  second  mode  has  a  characteristic  axial 
wavelength  which  is  twice  that  of  the  axisymmetric  mode.  Qualitatively  again,  this 
second  buckling  mode  is  similar  to  the  one  by  which  the  shells  failed  catastrophically 
in  the  bending  experiments. 

b.  26  <  D/t  <  40^ 

In  this  category  of  moderately  thick  shells,  the  prevalent  mode  of  instability  was,  again, 
short  wavelength  rippling.  Unlike  the  previous  category  of  shells,  a  well  defined  limit 
moment  was  recorded  in  the  moment-curvature  response  of  these  shells.  In  addition, 
some  non-uniform  growth  of  ovalization  was  observed  just  prior  to  failure  for  all  test 
specimens  in  this  category.  The  localized  deformation  was  typically  8-10  diameters 
^Approximate  values. 
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long.  In  a  few  cases,  more  than  one  area  of  localized  deformation  was  observed  to  grow 
simultaneously  (probably  triggered  by  initial  imperfections).  Pockets  of  short  wave 
length  ripples  were  usually  observed  in  the  trough  of  the  localized  region.  Following 
the  limit  load,  the  shell  buckled  catastrophically  by  the  appearance  of  one  sharp,  local 
kink  initiated  from  the  ripples. 

Two  possible  explanations  can  be  given  for  this  behavior.  The  first  is  based  on 
the  observation  that  the  curvature  at  rippling  was  progressively  closer  to  the  limit  load 
induced  by  uniform  ovafization.  Limit  loads  of  this  type,  which  we  will  call  natural,  are 
generally  not  very  imperfection  sensitive.  However,  as  will  be  shown  in  the  second  part 
of  this  work,  some  amplification  of  long  wave  length  imperfections  is  possible  as  the  limit 
load  is  approached.  Thus,  conceivably,  the  growth  of  long  wavelength  imperfections 
observed  in  the  experiments  was  due  to  this.  The  localized  ovalization  in  turn  triggered 
the  short  wavelength  ripples,  which  made  the  shell  axially  more  compliant,  and  caused 
the  limit  load  to  develop  earlier. 

A  somewhat  different  sequence  of  events  also  seems  possible  at  this  juncture.  It 
can  be  argued  that  the  instability  which  triggers  the  observed  behavior  is  still  rippling. 
For  this  range  of  Dji  values,  the  rippled  shells  remain  stable  longer  (i.e.,  the  difference 
between  and  is  larger)  than  thinner  ones.  The  ripples  grow,  make  the  shell  axially 
more  compliant,  and  cause  a  reduction  of  the  overall  bending  rigidity.  This  results  in 
the  development  of  the  limit  load.  The  presence  of  the  limit  load  is  then  responsible  for 
the  non-uniform  deformation  of  the  shell.  Some  support  for  this  argument  also  comes 
from  Calladine  (1982).  Although  he  did  not  address  localization  of  ovalization,  he  did 
in  essence  propose  that  the  ripples  cause  a  reduction  in  the  stiffness  of  the  shell  and 
the  development  of  a  limit  moment. 

These  questions  will  be  further  addressed  in  the  second  part  of  this  study  with 
the  help  of  numerical  results. 
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c.  D/t  <  26^ 


The  response  of  this  category  of  relatively  thick  shells  to  bending,  is  governed  by  the 
natural  hmit  load  instabihty  caused  by  uniform  ovalization.  Close  to  the  limit  load, 
long  wavelength  imperfections  tend  to  get  amplified.  The  ovalization  starts  to  grow 
non-uniformly  along  the  length  and  localizes  significantly  following  the  limit  load.  The 
shell  eventually  fails  in  the  region  of  localization.  A  preliminary  conclusion  is  that  the 
localization  is  a  direct  consequence  of  the  limit  load  instability  and  affects  the  behavior 
of  the  shell  primarily  after  the  limit  load.  Thus,  the  limit  load  calculated  by  assuming 
the  shell  to  ovalize  uniformly  is  very  representative  of  the  actual  value. 

The  bounding  Dji  values  of  the  three  categories  of  shell  behavior  are  only  ap¬ 
proximate  and,  very  definitely,  vary  with  the  shell  material  properties. 

Design  Considerations 

The  discussion  above  clearly  demonstrates  the  complexity  of  plastic  buckling  and  vari¬ 
ous  mechanisms  of  collapse  of  shell  structures.  Simulation  and  prediction  of  the  insta¬ 
bilities  discussed  is  of  great  fundamental  interest,  and  if  successfully  carried  out  should 
improve  the  general  understanding  of  the  problem  area.  However,  the  predictions  will 
obviously  be  complicated  by  the  three-dimensional  nature  of  the  mechanical  phenomena 
involved. 

It  is  possible  that  a  simpler  approach  may  adequately  serve  the  needs  of  structural 
design.  The  two  main  instabilities  observed  experimentally  are  the  short  wavelength 
rippling  and  the  limit  moment  resulting  from  uniform  ovalization.  Thus,  as  a  first  step, 
the  uniform  tube  analysis  of  Shaw  and  Kyriakides  (1985)  was  used  to  calculate  the 
limit  moment,  the  corresponding  curvature  and  Saint- Venant  bending  strain  for  shells 
with  D/ts  ranging  from  15  to  100.  The  average  values  of  the  key  material  parameters 
of  the  shells  tested,  and  listed  below,  were  used  in  these  calculations. 

E  =  10^  ksi(68.95  GPa),  Cy  =  42.6  ksi(293.7  MPa),  n  =  29  and  v  =  0.32  . 

Following  Reddy  (1979),  the  calculated  maximum  Saint-Venant  bending  strain 
^Approximate  value. 
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corresponding  to  the  limit  moment  is  plotted  against  Djt  (log-log  scales)  in  Fig.  24. 
The  critical  strains  measured  in  the  eleven  experiments  (i.e.,  either  ei  or  £;,)  are  also 
included  in  the  same  figure.  (Note  that  due  to  the  ovalization  of  the  cross  section 
the  actual  critical  axial  bending  strain  in  the  shells  is  between  2-3%  lower  than  the 
maximum  Saint- Venant  bending  strain.) 

An  independent  bifurcation  check  analysis  developed  by  Ju  and  Kyriakides  (1990) 
was  incorporated  in  the  uniform  tube  bending  einalysis  mentioned  above.  This  was  used 
to  identify  the  onset  of  ripple- like  bifurcations  during  the  bending  history.  The  buckling 


mode  used  is  as  follows: 


=  R  cos  px^^Cn  cos  nO  , 


V  =  cos  pa:  ^  sin  n0  ,  P=T —  ■<  (5) 

n=l 

Nu 

u  =  Rs'in  px  ^  e„  cos  n0  . 

n=0 

where  {u,v,  ti}  are  displacements  in  the  axial,  circumferential  and  radial  directions, 
measured  from  a  circular  toroidal  reference  shell  with  the  following  principal  curvatures 

KCOSO  1 

=  T1 — D - a  '^2  =  D  • 

1  -t-  kR  cos  6  R 

K,  in  eq.  (6)  is  the  curvature  of  the  bent  tube  at  the  time  of  the  bifurcation  test. 
The  formulation  used  is  based  on  Sander’s  (1963)  nonlinear  shell  kinematics.  As  is 
customary,  the  deformation  theory  of  plasticity,  with  the  state  of  stress  calculated 
in  the  tube  analysis,  was  used  in  the  bifurcation  check  together  with  Hill’s  concept  of 
comparison  solid  (see  Hill,  1958). 

The  critical  bifurcation  strains,  calculated  for  the  same  group  of  shells  mentioned 
above,  are  included  in  Fig.  24.  For  this  material,  the  calculated  bifurcation  strain  is 
lower  than  the  strain  corresponding  to  the  limit  load  for  shells  with  Djt  >  22.  For  shells 
with  lower  Dfl  values,  the  calculated  limit  load  precedes  the  bifurcation  instability.  The 
correlation  of  the  predictions  with  the  experimental  results  is  exceptionally  good.  A 
small  difference  in  the  value  of  the  transitional  Dfl  exists,  in  that  no  ripples  were 
observed  in  the  experiments  for  shells  with  Dft  <  25.7.  This  difference  is  partly  due 
to  variations  in  material  properties  from  shell  to  shell.  As  in  other  plastic  bifurcation 
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problems  (see  Hutchinson  (1974)),  the  corresponding  bifurcation  strains  calculated  by 
using  the  J2  flow  theory  are  unreasonably  high  for  thinner  shells  and  do  not  exist  for 
thicker  shells  in  the  range  of  D/ts  studied. 

The  same  degree  of  agreement  between  experimental  and  predicted  critical  strains 
htis  been  obtained  for  steel  shells.  However,  the  transitional  D/t  was,  in  general, 
higher  than  in  the  aluminum  shells.  Similar  success  of  this  simplified  approach  was 
also  demonstrated  for  the  problem  of  bending  under  external  pressure  in  Corona  and 
Kyriakides  (1988)  and  Ju  and  Kyriakides  (1990).  It  can  thus  be  concluded  that,  in 
spite  of  the  simplicity  of  the  approaches  followed,  the  limit  load  and  ripple  bifurcation 
predictions  constitute  very  useful  tools  for  establishing  the  limit  to  which  long  shells 
can  be  safely  bent. 

Another  simplification,  often  made  in  the  literature,  is  to  use  the  critical  strain  at 
which  axially-loaded  circular  cylindrical  shells  first  buckle  into  an  axisymmetric  mode 
as  the  critical  design  strain  of  shells  in  bending  (e.g.  Seide  and  Weingarden  (1961)  for 
elastic  shells;  Reddy  (1979),  Yun  and  Kyriakides  (1990)  for  elastic-plastic  shells).  The 
critical  axial  stress  of  this  problem  is  given  by 


<Te  = 


C'llC'22  —  Cv 


12 


1/2 


\d)  ■ 


The  corresponding  critical  axial  half  wavelength  is  given  by 


Ac  =  TT 


12(C„C22-C?2)J 


2  J 


{-) 


(8) 


where  Caff  are  the  instantaneous  moduli  of  the  material  at  bifurcation  (in  our  case, 
calculated  by  Jj  deformation  theory).  The  critical  strain  can  be  evaluated  from  (7) 
and  (3). 

The  critical  strains  predicted,  using  this  approach,  are  also  shown  in  Fig.  24.  As 
might  be  expected,  in  the  range  of  D/t  values  considered,  these  predictions  represent 
a  lower  bound  to  the  experimental  results.  For  lower  D/ts  the  predictions  are  low  by  a 
factor  of  two.  At  higher  D/t  values,  they  approach  the  bifurcation  predictions  obtained 
from  the  tube  analysis  mentioned  above.  In  view  of  this,  for  shells  with  D/t  <  80,  say, 
the  analysis  which  includes  the  effect  of  ovalization  should  be  prefered  over  eqs  (7)  and 
(8). 
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In  a  number  of  recent  publications,  investigators  have  explored  design  schemes  for 
assessing  the  post  bucking  strength  of  tubular  structures.  We  stress  that  pure  bending 
experiments  on  both  steel  and  aluminum  tubes  have  clearly  shown  that,  following  local 
buckling,  the  moment  and  deformation  capacity  of  shells  is  structurally  insignificant. 
Persistent  rotation  of  the  ends  of  the  shell  results  in  continued  drop  in  moment  and 
eventually  leads  to  fracture  and  total  failure.  This  is  due  to  the  very  localized  nature  of 
the  deformations  induced  to  the  shell  by  buckling.  Although  the  post  buckling  strength 
of  tubular  components  which  are  more  constrained  by  the  remainder  of  the  structure 
should  be  somewhat  better,  the  essence  of  the  statements  above  is  unaltered.  Thus  the 
post  buckling  regime  should  be  tackled  with  great  caution. 
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NOMENCLATURE 


shell  outside  diameter 

maximum  and  minimum  shell  diameter 

shell  mean  diameter 

Young’s  modulus 

shell  half  length 

moment 

aoDlt 

Ramberg- Osgood  hardening  parameter 
anisotropy  parameter 
shell  wall  thickness 

maximum  and  minimum  wall  thickness 


initial  ovality 


D  -  D  ' 


V  D^^  +  D^J 

£b,  £Li  £c  maximum  Saint- Venant  bending  strain  at  bifurcation,  limit  load  and 
collapse 

K  curvature 

Ki  tjDl 

«c  bifurcation,  limit  load  and  collapse  curvature 

A  axial  ripple  half-wave  length 

Ad  half- wave  length  of  predicted  ripples  (J2  def.  Theory) 

<7y  Ramberg-Osgood  yield  parameter 

cTo  yield  stress 


thickness  variation  parameter 


^max  ^min 
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Table  I  Geometric  and  material  characteristics  of  shells  tested 


Moment-curvature  responses  of  aluminum  tube 
("a"  s  limit  moment;  "i"  s  short  wavelength  rippling) 


Ovalization  along  shell  length  for  different  values  of  curvature 
(Aluminum-6061 -T6;  D/t  =  35.7.  Experiment  by  S.  Kyriakides 
Juarez,  1982,  unpublished) 


If.... 


Fig.  6  Stationary  ovaiization  transducer 


Pie.  7 


(a)  Transducer  for  scanning  ovalization  along  shell  length 

(b)  Transducer  schematic  and  A/D  convener  triggering  circuit 
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Uniform  Tube  Prediction 
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Bending  response  of  aluminum  shell  with  D/t  =  50 

(a)  Moment-curvature 

(b)  Ovalization-curvature 


o  > 
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Bending  response  of  aluminum  shell  with  D/t  =  32.2 

(a)  Moment-curvature 

(b)  Ovalization-curvature 


AI-606I-T6 


^=25.3 


Experiment 

Uniform  Tube  Prediction 


Bending  response  of  aluminum  shell  with  D/t  =  25.3 

(a)  Moment-curvature 

(b)  Ovalization-curvature 
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Bending  response  of  aluminum  shell  with  D/t  =  19.5 

(a)  Moment-curvature 

(b)  Ovalization-curvature 
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Fig.  24  Critical  maximum  St.  Venant  bending  strain  as  a  function  of  shell  D/t. 
Comparison  of  experimental  values  and  various  predictions 


